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We propose a prescription to quantize classical monomials in terms of symmetric and ordered 
expansions of non-commuting operators of a bosonic theory. As a direct application of such quanti- 
zation rules, we quantize a classically time evolved function 0{q,p,t), and calculate its expectation 
value in coherent states. The result can be expressed in terms of the application of a classical 
operator which performs a Gaussian smoothing of the original function O evaluated at the center of 
the coherent state. This scheme produces a natural semi-classical expansion for the quantum expec- 
tation values at a short time scale. Moreover, since the classical Liouville evolution of a Gaussian 
probability density gives the same form for the classical statistical mean value, we can calculate the 
first order correction in h entirely from the associated classical time evolved function. This allows 
us to write a general expression for the Ehrenfest time in terms of the departure of the centroid of 
the quantum distribution from the classical trajectory provided we start with an initially coherent 
state for each subsystem. In order to illustrate this approach, we have calculated analytically 
the Ehrenfest time of a model with coupled non-linear oscillators with non-linearity of even order. 

Pacs numbers: 03.65.Sq, 05.45.-a, 05.20.Gg 
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I. INTRODUCTION 

Since the earliest days of the quantum theory, the in- 
vestigation of the differences between the probabilistic 
concepts in classical Liouville and quantum dynamics has 
been an important issue. There has been many stud- 
ies in the last two decades concerning the semi-classical 
regimen, of systems whose classical counterparts exhibits 
chaos l] . The question of estimating how long the classi- 
cal and quantum evolution stay close has been one of the 
main questions in semi-classical analysis. For classically 
chaotic flows, there was a conjecture which is now rig- 
orously proved "s] , that the break time or Ehrenfest time 
{ts) diverges logarithmically with h. In the classically 
regular flow, it was suggested in ,4] that the behavior of 
ts with h is algebraic {h~^ ), but no universal nature of 
such behavior has been shown yet. 

According to the famous Ehrenfest 's work [^, for quan- 
tum states which are localized enough, the time variation 
of the mean quantum momentum must be equal to the 
local force. This statement is exact for quadratic Hamil- 
tonians, but its validity is restricted to a short time scale, 
the Ehrenfest scale, for more general non-linear systems. 
Mathematically, the Ehrenfest theorem allows us to write 
{0{X)) = 0{{X)) for times smaller than the Ehrenfest 
time. In this situation, the initial dynamics is described 
essentially by a mean field approximation, where we have 
a localized packet obeying classical equations of motion. 
Then, it is reasonable to expect a rather good agreement 
between quantum and classical Liouvillian centroids and 
the classical trajectory. In fact, this scenario have already 
been reported in the literature 0. 

In this work, we propose a simple analytical scheme to 
calculate the Ehrenfest time for integrable systems. Our 
starting point is to propose a classical Liouvillian opera- 
tor that makes explicit the symmetric form of the usual 



quantization rules. Using such an operator, we are able to 
take a general classical function O which expresses the 
time evolution of a physical quantity, quantize it, and 
evaluate its expectation value in coherent states at each 
time. The final result, which is shown to be analytically 
identical to the statistical average calculated through the 
classical Liouville formalism, is written in terms of a cer- 
tain differential operator acting onto the original classical 
function. Since this is just a compact form of express- 
ing the action of the corresponding power series in h, 
our recipe automatically leads to a semi-classical expan- 
sion around the time evolved classical function. Notice 
that the equality between quantum and statistical cen- 
troids guarantees that we are working within a classical 
time scale. This allows us to define mathematically the 
Ehrenfest scale using just the Liouvillian Gaussian aver- 
age without solving the quantum dynamical problem. 



As an example, we present an explicit calculation of the 
Ehrenfest time for iV-coupled non-linear oscillators, with 
nonlinearity of order 2fc, an even integer number. This 
model to which we will address the quantum-classical 
departure issue can be associated to several nonlinearly 
interacting fields via Kerr- type 0,01 and cross-Kerr type 
interaction fl(| known to be relevant in quantum optics 
and also for two quantized vibrational modes of a 
single trapped ion It is an integrable model where 
the role of nonlinearity can be studied analytically for 
several quantities, and, in particular, for N = k = 2 the 
quantum-classical break time has been determined [T^ 
based on physical properties of the exact quantum states. 
Although the emphasis of studying the break times in the 
literature has been on 'chaotic states' 0,01, iieie we are 
concerned with integrable cases, where we are able to 
derive an analytical expression for the Ehrenfest time. 
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II. ORDERED QUANTIZATION 



B. Normal Ordering in Bosonic Operators 



A. Definition 

Let us start by presenting a convenient quantization 
scheme for a single degree of freedom system, which wih 
be easily generalized to degree N. Consider two non- 
commuting operators A and B whose commutator is a 
c-number, denoted by c ([A, B] = c). Then, a given 
classical monomial a"6™, where a and b are canonically 
conjugated classical variables ^2^, will be quantized in a 
symmetric and ordered form, through the prescriptions: 



a^b"" S^ g , (1) 



where the super-operator Sx „ is given by 



S 



A,B 



k=0 



fc! 



(2) 



The above index (A, B) corresponds to the ordering with 
A on the left and B on the right. Since [9^,9^] = 

[9^,-8] = [A, 9^] — Q, S^j^ and Sj^ ^ are classical uni- 
tary differential super-operators, satisfying 5^ gSj^ ^ = 
Sg B = 1- The prescription defined in expression 
leads to different orderings (depending on which variable 
is chosen to be a or b) for the same original classical func- 
tion and, therefore, the associated symmetric operators 
must be the same. From this consideration, we deduce 
the following ordering formulas: 



B.A 



B"^A^ = S^- -A^B"^ = p'^^'^^^^Sb A^B^' 

A,B 



(3) 



As an immediate application of the formulas presented 
above, we will express the totally symmetric ordered ex- 
pression Sq pQ^P"^ in a form more suitable for our pur- 
pose of takmg expectation values in the Weyl-Heisenberg 
coherent states. We start by expressing such a polyno- 
mial in terms of bosonic creation and annihilation oper- 
ators 

SqpQ^P"' = e~'*^o9pg"p™ (4) 

n+m 



where Znm 
that 



(\ 'n-\-'rn 
y/hpA . Now, using the fact 



(a±at)" = ^(±ir^(")e^^^^at (fit ,(5) 
we can re-express (0J as 

n m 

with the sub-indexes that we introduced, to indicate 
where the action of the differentiation should take place. 
At the end of calculation we must erase all these sub- 
indexes. Using the relation ©, we get 



\rn — l , 



X 



k=0 1=0 

st 3^443 t 6 
°ie 3e 



(7) 



at f-in-k-i 

-3 l^a[ ai^ 



fA"^-l-f,l f.k 



This is the normal ordered expression (in the creation 
and annihilation operators) for the original QP-ordered 
monomial. Now, it is a simple matter to calculate its 
expectation value. 



Using a classical displacement operator, e°^="/(a;) = 
f{x + a), which can be used to write e"^^-*^^ ^"5™ = 
{A — ^dp)"B"\ one can show that these results repro- 
duce those in the textbooks (see e.g. Louisell ^3)- ^ 
simple example shows that our quantization scheme leads 
to the usual quantization rules, but in an automatically 
ordered form. Consider the following product of classi- 
cal canonical phase space variables: q^p. According to 
the usual ru les, this is transformed in a totally symmet- 
ric operator ^ of the form ^{Q^P + PQ^ + QPQ); and 
using the commutation relation [Q, P] = ih, re-written 
as Q^P — thQ. But this is exactly the result produced by 
the expression in ^ with the choice (a, 6) — {q,p), just 
by making some few derivatives. We finally note that the 
super-operators p has already appeared in the liter- 
ature in the case oi the canonical phase space operators 
{A, B) = (Q, P), in slightly different contexts 



C. Coiierent States Representation 

Since we are interested in the connection between 
quantum and classical mechanics, the coherent state ba- 
sis appears as the most appropriated one. In particular, it 
will be of particular interest for us to evaluate the expec- 
tation value in coherent states of some operator products 
like the one treated in previous section. Then, we first 
calculate the matrix elements in the coherent state basis 
of the operator function given in Eq. Q 



(ai Q!2) 



X e2 



dadd 



(a + 6) 



h{c-d) 



.(8) 
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Setting now ai = a2 = ao = and performing 

adequate variable transformations, we finally get 



ao|5QpQ"P"|ao> 



where 



i<lo,Po) = (^{aa\Q\ao) , {(^o\P\aa] 



(9) 



(10) 



The expectation value given in must be calculated 
through the series expansion of the classical operator 
6*^0, which gives a natural expansion in powers of h, 
showing its semi-classical nature. In fact, in the classi- 
cal limit h —f 0, the quantum expectation value of the 
operator function reduces to a purely classical function 
calculated at the center of the coherent packet located at 
(go,Po)- 

Another interesting relation can be obtained from a 
similar calculation 



(11) 



which implies that 



= ei^o ( ""H^"''''" ) ggpg^. (12) 

It is important to note that the term within the paren- 
theses in Ea. (|ll(l is exactly the Weyl transform of the 
operator The extra operator factor e9^° is 

what makes reference to the width of the coherent packet, 
as we shall see latter. These results also point out to the 
existence of an asymptotic classical limit {h — s- 0) for such 
expectation values in coherent states. 



III. SHORT TIME QUANTIZATION 

The usual quantization rules are defined in the Heisen- 
berg picture, where the Hamilton's equations solution 
q-Hiq,P,t) and pT-i{q,PTt), are transformed into Heisen- 
berg operators Qh{Q, t) and Ph{Q, P, t), where Q and 
P denote Schrodinger operators. On the other hand, 
since the Heisenberg and Schrodinger pictures coincide 
at t = 0, any scheme of quantization based on the 
Schrodinger picture would describe reasonably the quan- 
tum world for very short times. However, as long as we 
are interested in analyzing the quantum operator evolu- 
tion only during a classical (mean field) time scale, this 
would suffice. 

In this context, consider a classically time evolved func- 
tion O that has the following well defined expansion 



Oiq,p,t) 



n.m— 



(i)g>" 



(13) 



The coefficients Cn,m{t) contain all the time dependence. 
Now, applying the operator of ordered quantization Sq p 
and taking the average in coherent states we obtain 

{ao\d{Q,P,t)\ao) ^ei^oO{qo,po,t), (14) 

where we have used the result In order to estimate 
for how long the result f(T^ could be trusted, we present 
the calculation of the classical statistical counterpart of 
the problem. In the classical Liouvillian formalism, the 
following mean value is defined in the phase space (g,p) 

(O(go,po,t)) - J dqdp piq,p,t)Oiq,p,0), (15) 

where (goiPo) stands for the center of the initial distri- 
bution p{q,p,0). Then, since p{q,p,t) = e^*^ p{q,p,0) and 
0{q,p,t) — e^^*O{q,p,0), taking into account the fact 
that the volume in the phase space is preserved, we can 
re-write the previous equation in a Heisenberg-like form 

{Oiqo,po,t)) = j dqdpp{q,p,0)O{q,p,t). (16) 



Consider now a Gaussian initial distribution with width 
cr, and centered at the point (go,Po)- Performing the 
variable transformations, {q — qo) — x and (p — pq) = y, 
and using the classical displacement operator again, we 
can rewrite Ea. (|16|l as 



(O) = / dxdy 



3^^'e''^''C'(g,p,t). (17) 



Notice that the function O inside the integral no longer 
depends on the integration variables. Hence, by complet- 
ing squares and performing formally the integration, we 
finally obtain 



{O{qo,po,t)) =eTVoO(qo,po,t). 



(18) 



Now, we have proved that the effect of applying the op- 
erator eT^o is exactly that of smoothing the function O 
through a Gaussian mean, where a is related to the width 
of the Gaussian distribution to be used in the smooth- 
ing. Accordingly, this operator has already been used to 
express the Q-function as the Gaussian smoothing of the 
Wigner function It is important to notice that this 
result is exact for Gaussian statistical averages, i.e., it 
was derived from first principles without any approxima- 
tion. The only implicit assumption used in the calcula- 
tion was the existence of the derivatives in all orders for 
the function O. 

The comparison between expressions (|14() and H18|l . 
with a — h, confirms the fact that our proposal of a 
Schrodinger quantization for classical function should be 
adequate either at a short time scale or in cases in which 
the classical function depends linearly on the phase space 
coordinates {qo,po) (e.g. harmonic oscillator). The re- 
sults ()14|l and (|18|) can easily be extended for higher de- 
grees of freedom. 
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At this point one might formulate a question. First, 
we interpret the action of the operator e~^o on any clas- 
sical function, as an exact factorization of the effect of 
a Gaussian wave packet contribution. In fact, for com- 
pletely localized distribution {h — 0) there are no correc- 
tions. Then, one may ask: is it possible, in the partic- 
ular case of coherent separable initial states, to express 
the exact quantum centroid, (R.(t)) = in 

terms of the smoothed form e^^^rdt), where rdt) = 
{qc{t)TPc{ty) is a certain classical dynamical vector in 
phase space? The answer will be positive if we are able 
to calculate the inverse of the Gaussian smoothing oper- 
ation on the quantum centroid vector in the phase space 



rc(i) 



ao|R(i)|ao) 



(19) 



This function Vcit) will then be a coherent quantum tra- 
jectory, in the sense that it will carry all quantum dynam- 
ical information possible to put in a trajectory, except the 
contribution due to a Gaussian smoothing effect. For the 
trivial case of N non-interacting harmonic oscillators and 
the case of a bilinear coupling between two harmonic os- 
cillators, the quantum coherent and classical trajectories 
coincide ( (qdt) , pc{t)) = {qo{t),Po{t)) ), as expected. 



IV. THE EHRENFEST TIME 

A. Formal Definition 

Now, we have all the necessary tools to undertake the 
problem of estimating the Ehrcnfcst time in the case of 
an initially separable coherent Gaussian wave packet. As- 
suming that the break instant occurs when the first order 
correction in h becomes as important as the original vec- 
tor, we apply the smoothing process to the phase space 
classical trajectory vector of a system with N degrees of 
freedom r(i) = {qi{t),pi{t), - ■ ■ ,qN{t),pN{t)). Expand- 
ing the smoothing operator up to the first order in h, we 
obtain 



r{t) 



r(t) + -V\(i), 



(20) 



where 



= X^ili ^1 the A^-dimensional Laplacian 



operator. Now, we formally define the Ehrenfest time 
is as being the instant at which the magnitude of the 
difference between the quantum centroid and the corre- 
sponding classical vector in phase space becomes equal 
to the magnitude of the initial classical vector. Mathe- 
matically this condition can be expressed as follows: 



(k{tE)) - r{tE) 



= 1. 



(21) 



Finally, using Eq. H20() , we obtain the prescription for the 
analytical calculation of the Ehrenfest time 



h\\V^r{tE) 
4 1 



(22) 



It is remarkable that we do not need to solve the quantum 
equations of motion to find tE- In what follows, we will 
calculate the above defined break time for the system of 
N coupled non-linear oscillators. 



B. Application to a Non-Linear System 



Consider the following classical Hamiltonian 



N 



■ N 

E 

.1=1 



PI + Ti 



(23) 



where fc > 1 is integer and g is the only coupling constant 
of the system, and from which we define the characteris- 
tic classical action A = ^ (j^^~T^^ 01 • "^^^ equations 

of motion can be solved by noticing that A itself is a 
constant of motion. The result reads 




FIG. 1: Departure between quantum and classical centroids 
as a function of the dimensionless parameter cut for the bi- 
dimensional quartic oscillator {N = 2, fc = 2). These calcula- 
tions were analytically performed with qi = pi = 1 {A = 2), 
LJi = LJ = 1, g = 0.1 and h = 1 (solid line), h = 0.1 (dashed 
line) and h = 0.01 (dot-dashed line). The Ehrenfest time scale 
12511 for each value of h is represented by a vertical line in the 
corresponding style. 



(24) 



where {qi,Pi) are the initial conditions of the i-th oscil- 
lator and Qiit) — uoit + gtkK^^^ is a rotation angle in 
phase space. Finally, using 122|) and l|24|l . we get for the 
Ehrenfest time of this system 







COS Qi{t) 


sin9i(t) 










_ -sine.,(t) 


cos Qi{t) 







tE = 



1 



k{k-l) 



2A 



0.5 



hp' 

8A 



,(25) 
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where we kept only the first order terms in h. The cases 
(A^ = l,k = 2) and {N = 2,k = 2) can be seen to re- 
produce respectively the results calculated in ref. |23| and 
[2]|. Trivial limits are also contemplated by the above 
result, namely, the cases of an harmonic system {g = 
or fc = 1), for which ts oo. Moreover, our result is in 
accordance with some conjectures predicting the general 

form (t) for the break time of classically integrable 
systems [23, where fl = fc(fc_i)"gA<.-i and S = 2A are 
respectively the typical frequency and classical action of 
the system in consideration. We illustrate in figurcnthe 
Ehrenfest scale predicted by the expression H25|l for the 
case {N = 2, k = 2), where we can see an algebraic t^) 
short time departure. We also note that the first order 
correction to the Ehrenfest time scale seems to indicate a 
more appropriated parameter to measure the classicallity 
in this system: ^ <C 1. 

V. CONCLUSION 

In this paper we proposed a scheme for an ordered sym- 
metric quantization based on the action of certain differ- 
ential operators. As an application, wc quantized a clas- 



sically time evolved dynamical functions of the canonical 
variables in phase space, and showed that such a pro- 
cedure is adequate during the Ehrenfest time scale for 
separable coherent initial states. This allowed us to pro- 
pose a formal definition for the Ehrenfest time in terms 
of the phase space Laplacian operator acting on the clas- 
sical solutions of the equations of motion. We performed 
an explicit calculation for a system of N coupled nonlin- 
ear (2A:-th order) oscillators and calculated the Ehrenfest 
time for general N > 1 and A: > 1. The results are shown 
to agree with the results known in the literature for some 
particular cases. 
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